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Properties of the intracluster medium of Abell 3667 observed with Suzaku XIS* 
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Abstract 

We observed the northwest region of the cluster of galaxies A3667 with the Suzaku XIS instrument. The temper- 
ature and surface brightness of the intracluster medium were measured up to the virial radius (raoo ~ 2.3 Mpc). The 
radial temperature profile is flatter than the average profile for other clusters until the radius reaches the northwest 
radio relic. The temperature drops sharply from 5 keV to about 2 keV at the northwest radio relic region. The 
sharp changes of the temperature can be interpreted as a shock with a Mach number of about 2.5. The entropy slope 
becomes flatter in the outer region and negative around the radio relic region. In this region, the relaxation timescale 
of electron-ion Coulomb collisions is longer than the time elapsed after the shock heating and the plasma may be out 
of equilibrium. Using differential emission measure (DEM) models, we also confirm the multi-temperature structure 
around the radio relic region, characterized by two peaks at 1 keV and 4 keV. These features suggest that the gas is 
heated by a shock propagating from the center to outer region. 

Key words: galaxies: clusters: individual (Abell 3667) — X-rays:ICM, shock wave 
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1. Introduction 

Clusters of galaxies grow through gravitational infall and 
mergers of smaller groups and clusters. Merger events con- 
vert kinetic and turbulent energy of the gas in colliding sub- 
clusters into thermal energy by driving shocks in the cluster. 
A fundamental issue of cluster growth is the nature of heating 
processes in merger events. Although shocks caused by merger 
events are very important for gas heating and particle acceler- 
ation, there are only a few clusters for which clear evidence 
of shocks has been obtained (1E0657-56, A520: Clowe et al. 
2006; Markevitch et al. 2005). 

An additional science issue of merging clusters is the tem- 
perature structure. In dynamically young systems, the intra- 
cluster medium (ICM) is thought to be in an early stage of ther- 
mal relaxation. Some theoretical studies predict non-ionization 
equilibrium states and an electron-ion two-temperature struc- 
ture in the ICM of merging galaxy clusters (Rudd & Nagai 
2009; Akahori & Yoshikawa 2010). Recent Suzaku results 
show that the temperature profile of the merging cluster A2142 
is slightly in excess compared to other relaxed clusters and a 
possible deviation between electron and ion temperatures may 
be present (Akamatsu et al. 2011). However, there is little 
information about the ICM properties in the outer regions of 
merging clusters because of the faintness of the X-ray surface 
brightness. 

In this paper, we present results from observations of the 
merging cluster A3667 (z = 0.0556) with XIS on board Suzaku 
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(Mitsuda et al. 2007). It is a very bright merging cluster with 
irregular morphology and an average temperature of 7.0 + 0.5 
keV (Knopp et al. 1996). Briel et al. (2004) studied the ICM 
characteristics in the cluster central region and showed that the 
ICM emission is elongated to the northwest direction. The 
most striking features of A3667 is a discontinuity in the X- 
ray surface brightness, namely a "cold front", which was found 
with Chandra (Vikhlinin et al. 2001). Another feature of A3667 
are the two extended, symmetrically located regions of dif- 
fuse radio emission revealed with 843 MHz ATCA observatory 
(Australia Telescope Compact Array: Rotttgering et al. 1997). 
The presence of these structures naturally indicates that A3667 
is a merging system and the subcluster infall seems to be oc- 
curring along the northwest-southeast direction. Actually, the 
overall cluster emission is elongated in this direction. 

A recent XMM-Newton observation of the northwest radio 
relic of A3667 has shown significant jumps in temperature and 
surface brightness across the relic (Finoguenov et al. 2010). 
The observed ICM properties associated with the jump indicate 
the existence of a shock front with a Mach number Ai ~ 2.4. 

From Suzaku XIS/HXD observations of A3667, Nakazawa 
et al. (2009) measured the cluster emission up to 40' and set 
upper limits on the non-thermal X-ray emission in XIS spectra 
of 7.3 X 10 erg cm"^ s ' when extrapolated to 10-40 keV. 
Nakazawa et al. (2009) also observed very high temperature 
components {kT ~ 14 keV) in the cluster center Since their 
main interest was to look into the non-thermal phenomenon, 
the global ICM properties in the cluster outskirts remain to be 
further explored. 

The purpose of this study is to clarify the physical state of 
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Table 1. Suzaku observation log of Abell 3667 



Region(Obs.ID) 


Obs. date 


Raw Exposure* 


Cleaned Exposure ' 






(ks) 


(ks) 


CENTER (801055010) 


2006-5-6T17:40:17 


19.3 


16.8 


OFFSET 1 (802030010) 


2006-5-6T07:03:14 


16.2 


11.8 


OFFSET2 (802031010) 


2006-5-3T17:47:01 


81.3 


60.3 



*:C0R2 > GV t:C0R2 > 8 GV and ELV 10° 





0.5 



IliiiMt 



1 2 

Energy (keV) 



10 



Fig. 1. X-ray image of the northwest region of A3667 in the energy 
band 0.5-8.0 keV, after subtraction of the NXB with no vignetting cor- 
rection and smoothing by a 2-dimensional gaussian with o" = 1 6 pixel 
= 17". The large green circle indicates the virial radius of A3667, and 
black dotted circles show the annular regions used for the spectral anal- 
ysis. Small green circles show point sources that are detected with 
XMM-Newton. The SUMMS 842 MHz radio image is shown with 
magenta contours. 

the ICM in the outer regions of the merging cluster A3667. 
We use Ha = 70 km s"' Mpc-\ = 0.27 and Qa = 0.73, 
respectively, which means that 1' correspond to a diame- 
ter of 66 kpc. The virial radius is approximated by raoo = 
2.77/!7„i(<r)/10keV)'/2Mpc/£(z), where E{z) = (Qm(1 +z)^ + 
1 - i^ivi)'^^ (Henry et al. 2009). For our cosmology and red- 
shift, r2oo is 2.26 Mpc (= 34.' 1) with kT = 7.0 keV. In 
this paper, we employ solar abundances given by Anders & 
Grevesse (1989) and Galactic hydrogen column density of 



A^H = 4.7 X 10^" cm--' (Dickey & Lockman 1990). Unless oth- 
erwise stated, the errors correspond to 90% confidence for a 
single parameter 

2. Observations & Data reduction 

As shown in Fig 1, Suzaku carried out three pointing obser- 
vations of Abell 3667 along the northwest merger axis in May 
2006, designated as Center, Offsetl, and Offset2. The obser- 
vation log is summarized in Table 1. All observations were 
performed with either normal 5 x 5 or 3 x 3 clocking mode. 

The combined observed field extends to the virial radius of 
A3667 (34.1 ~ 2.26 Mpc). During the observations, the con- 
tamination from the Solar Wind Charge eXchange (SWCX) 
was not significant, which we base on ACE SWE data and 



Fig. 2. The spectrum of the outermost region used for the background 
estimation (r = 37.'8 - 42.'0), after NXB subtraction. The XIS BI 
(Black) and FI (Red) spectra are fitted with CXB + Galactic compo- 
nents (LHB, MWH) {apec+wabs(apec-¥powerlow}). The CXB spec- 
trum is shown with a black curve, and the LHB and MWH components 
are indicated by green and blue curves, respectively. 

the absence of flare-like features in the XIS light curve (see 
Appendix 1). 

The XIS instrument consists of 4 CCD chips: one back- 
illuminated (BI: XISl) and three front-ifluminated (FI: XISO, 
XIS2, XIS3) ones. The IR/UV blocking filters had a signif- 
icant contamination accumulated at the time of the observa- 
tions, and we included its effect and uncertainty on the soft 
X-ray effective area in our analysis. We used HEAsoft ver- 
sion 6.9 and CALDB 2010-12-06 for ah the Suzaku data anal- 
ysis presented here. In the XIS data analysis, we limited the 
Earth rim ELEVATION > 10° to avoid contamination of scat- 
tered solar X-rays from the day Earth limb. Additionally, we 
performed event screening with cut-off rigidity (COR) > 8 
GV to increase the signal to noise ratio. We extracted pulse- 
height spectra in 10 annular regions whose boundary radii were 
4.'2, 8.'4, 12.'6, 16.'8, 21.'0, 25.'2, 29.'4, 33.'6,37.'8 and42.'0, 
with the center at (20''12'"3P,-56°49"12^). We analyzed the 
spectra in the 0.5-10 keV range for the FI detectors and 0.5- 
8 keV for the BI detector. In all annuli, positions of the cali- 
bration sources were masked out using the calmask calibration 
database (CALDB) file. 

3. Background modeling 

An accurate measurement of the background components 
is important for the ICM study in the cluster outer region. 
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Table 2. Estimation of tlie CXB fluctuation. 



Region 


0'-4.'2 


4:2-8:4 


8:4-12:6 


12:6-16:8 


16:8-21:0 


21:0-25:2 


25:2-29:4 


29:4-33:6 


33:6-37:8 


SRR' 


22.28 


20.89 


3.53 


0.42 


0.29 


0.16 


0.05 


0.04 


0.03 


e,Suzaku 


50.3 


158.9 


73.2 


67.7 


80.3 


72.1 


59.9 


62.5 


60.0 


c/^cxbI 


14.0 


7.9 


11.6 


12.1 


11.1 


11.7 


12.8 


12.6 


12.8 



* : S ROUS JiATIO^EG[%]:The fraction of simulated cluster photons that full in the region compared to the total photons 
generated in the entire simulation. 

t:[arcmin^] 

$:[%]: In this paper, := 5 x 10"'"* erg cm"^ s^'is assumed for all regions. 5^ is the flux limit of the extracted point sources. 



Table 3. Best-fit background parameters 







NOMINAL 


CXBMAX 


CXBMIN 


CONTAMI+10% 


CONTAMI-10% 


LHB 














kTikeV) 




0.11 ±0.01 


0.11 ±0.01 


0.10 ±0.01 


0.11 ±0.01 


0.10 ±0.01 


norm* (xlO" 




4 85+1-™ 


4 48+1.51 


5 59+1-5'' 

-^•-^^-1.80 


4 25+1-30 


6 24+1-^1 


MWH 














kTikeV) 




0.37 ±0.03 


0.37 ±0.04 


0.37 ± 0.04 


0.37 ±0.04 


0.36 ±0.04 


norm* (xlO" 




11.6!!-^ 




11 3+'-2 


1 1 o+'i 




^7d.o.f 




381/344 


398/344 


406/344 


381/344 


382/344 



*:Normalization of the apec component scaled with a factor l/400w. 



Norm=5^ Jn^n/f JV/(4w(l + z^)D^) x 10 '^cm ^arcmin ^, where Da is the angular diameter distance to the source. 



We examined four background components, namely non-X-ray 
background (NXB), cosmic X-ray background (CXB) and the 
Galactic emission consisting of the Milky Way Halo (MWH) 
and the Local Hot Bubble (LHB). The NXB component 
was estimated from the dark Earth database by the xisnxbgen 
FTOOLS (Tawa et al. 2008) and was subtracted from the data 
before the spectral fit. To adjust for the long-term variation of 
the XIS background due to radiation damage, we accumulated 
the NXB data for the period between 150 days before till 150 
days after the observation of A3667. 

3.1. Cosmic X-ray Background Intensity and Fluctuation 

Using seven XMM-Newton data sets reported in Briel et al. 
(2004) and Finoguenov et al. (2010), we searched for point-like 
sources within the Suzaku field of view using the wavdetect 
tool in the CI AO package. As shown by the green circles in fig- 
ure 1 , we detected and subtracted 22 point sources which have 
2.0-10 keV fluxes higher flian 5 ^ = 5 x lO"^"* erg cm^^ s ' . The 
extraction radius was set to 1' or 2' corresponding to the HPD 
of the Suzaku XRT. We estimated the CXB surface brightness 
after the source subtraction to be 5.97 x 10"*^ erg cm"^ s"^ sr~i 
based on ASCA GIS measurements (Kushino et al. 2002). 

Kushino et al. (2002) set the flux-limit of eliminated point 
sources to be 5'c = 2x10"'^ ergcm"^ s we subtracted sources 
down to 5 X 10 erg cm"^ s ' . This causes our flux-limit value 
to be sufficiently lower than the Kushino et al. level. 

We need to estimate the ampUtude of the CXB fluctuation, 
caused by the statistical fluctuation of point sources in the FOV. 
We scaled the measured fluctuations from Ginga (Hayashida 
1989) to our flux limit and field of view foUowing Hoshino et 
al. (2010). The fluctuation width scales as -0.5 power of the 
effective field of view {Q.^) and 0.25 power of the point-source 
detection Umit. We show the CXB fluctuation level, and 



o'/IcxB, for each spatial region in table 2. 

3.2. Galactic Components 

To estimate the Galactic background components, we use 
the data in the outermost region which is located outside of 
the virial radius (37.'8 - 42:0). We account for the Local 
Hot Bubble (LHB: ~ 0.1 keV with no absorption), the Milky- 
Way Halo (MWH: ~ 0.3 keV with absorption) and the CXB 
component mentioned above. The model is described by 
apec + wabs(apec + powerlaw). The redshift and abundance in 
both the apec components were fixed at and unity, respec- 
tively. To represent the Galactic and CXB emission, we used 
the XIS response for a source of uniform brightness. The spec- 
trum was well fitted with the above model (xl^j^l.lQ for 344 
degrees of freedom). The temperatures of the LHB and the 
MWH are 0.11 ± 0.01 keV and 0.37 ± 0.03 keV, respectively, 
consistent with the typical Galactic emission. The resultant pa- 
rameters and the spectrum are shown in Table 3 and in Figure 2. 

3.3. Stray light 

Stray light consists of photons entering from the detec- 
tor outside the FOV. The Suzaku optics often show signif- 
icant effects of stray light from nearby bright X-ray sources 
(Serlemitsos et al. 2007). The extended point spread function 
of the telescope with a half -power diameter of 1 .'7 also causes 
contamination of photons from nearby sky regions. 

To examine the effect of stray Ught from the bright central 
region of the A3667 cluster, we calculated the photon contribu- 
tions to our annular regions from different regions of the clus- 
ter using the ray tracing code xissim (Ishisaki et al. 2007). We 
used version 2008-04-05 of the simulator. The surface bright- 
ness distribution is one of the input parameters necessary to run 
xissim. We used a yS-model (fi = 0.54, = 2.'97) based on the 
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Table 4. Relative count contributions (%) due to PSF broadening of the Suzaku XIS mirror. Relative contri- 

bution (%) from each sky region in the annular extraction region due to the PSF broadening of the Suzaku 
XIS min'or. Most of the PSF contaminations in an annular originates from the inner neighboring annular. 



Detector/ Sky 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(l)0'-4.'2 


90.9 


8.9 


0.1 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


(2)4.'2-8.'4 


20.6 


75.4 


3.9 


0.1 


0.0 


0.0 


0.0 


0.0 


0.0 


(3)8.'4 - 12.'6 


3.0 


21.0 


71.9 


3.7 


0.2 


0.1 


0.0 


0.0 


0.0 


(4)12.'6 - 16.'8 


0.9 


1.4 


12.7 


75.1 


9.3 


0.4 


0.1 


0.1 


0.0 


(5)16.'8-21.'0 


0.8 


0.6 


1.4 


14.3 


73.8 


8.4 


0.4 


0.1 


0.0 


(6)21.'0-25.'2 


1.1 


1.2 


0.7 


1.5 


14.1 


72.1 


8.6 


0.3 


0.2 


(7)25.'2-29.'4 


1.9 


1.4 


0.9 


0.8 


1.2 


13.3 


72.1 


7.7 


0.5 


(8)29.'4-33.'6 


0.8 


0.6 


0.5 


0.4 


0.4 


0.9 


13.4 


73.4 


9.1 


(9)33.'6-37.'8 


1.8 


1.6 


1.1 


0.7 


0.4 


0.7 


0.5 


11.2 


71.4 



(a)0'-4.'2 



(b)4.'2-8.'4 



(c)8:4-12.'6 




0.5 

(d)12.'6-16.'8 



Energy (keV) 




0.5 

(g)25.'2-29.'4 



Energy (keV) 




Energy (keV) 




0.5 

(e)16.'8-21.'0 



Energy (keV) 




0.5 

(h)29.'4-33.'6 



Energy (keV) 




Energy (keV) 



< csj 'T 





— 1 — 1 — 1 — 1 1 1 1 • 1 











0.5 

(f)21.'0-25.'2 



Energy (keV) 



0.5 

(i)33.'6-37.'8 



2 5 
Energy (keV) 



... '..::iJ 



Energy (keV) 




Fig. 3. NXB subtracted spectra in each annular region. The XIS Bl (Black) and Fl (Red) spectra are fitted with the ICM model {wahs + apec), along 
with the sum of the CXB and the Galactic emission {apec + wabs(apec + powerlaw}). The CXB component is shown with a black curve, and the LHB 
and MWH emissions are indicated by green and blue curves, respectively. The total background components are shown by the orange curve. 
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Table 5. Best-fit parameters of the ICM 





0'-4.'2 


4r2-8.'4 


8^4 - 12.'6 


12^6 - 16.'8 


16.'8-21.'0 


21.'0-25.'2 


25.'2 - 29.'4 


29.'4 - 33^6 


33.'6-37.'8 


NOMIMAL 


kT{keV) 


7 14+u.n 

'■^^-0.17 


6 57+" " 




6 94+" " 


6 49+" ''^ 


5 70+"-^'* 

■^■'"-0.66 


c -Ji+l.ll 

J.JO Q 


1 92+'"^ 
'■^^-0.37 




Z 


0.29 ± 0.03 


0.33 ± 0.03 


0.26 ±0.09 


0.20 ±0.11 


0.22 ±0.12 


0.26 ± 0.22 


0.2 (FIX) 


0.2 (FIX) 


0.2 (FIX) 


norm* 


631.2±6.3 


209.1 ± 1.8 


66.2 ± 1.9 


45.7 ± 2.0 


27.3 ± 1.4 


16.3 ± 1.4 


5.0±0.5 


2.9 ±0.4 


1.0±0.3 


cll 

'^0.4-lOkcV 


91.80±0.69 


29.54 ±0.19 


10. 13 ±0.66 


6.46 ± 0.06 


3.67 ±0.07 


2. 10 ±0.05 


0.63 ±0.01 


0.14 ±0.01 


0.04 ±0.01 


X^/d.o.f 


386/31 1 


39 1/311 


323 / 311 


498 / 404 


371 / 330 


348 / 330 


404 / 347 


388 / 347 


454 / 347 


CXB MAX+NXB 3'/< ADD 


kT{keW) 


7.11 


6.51 


6.76 


6.91 


6.22 


5.40 


4.10 


0.99 


1.05 


Z 


0.29 


0.33 


0.25 


0.19 


0.21 


0.24 


0.2 (FIX) 


0.2 (FIX) 


0.2 (FIX) 


norm 


612.3 


201.8 


63.6 


46.1 


26.0 


15.4 


4.4 


2.1 


0.7 


s" 


88.85 


28.47 


9.63 


6.45 


3.44 


1.92 


0.47 


0.03 


0.01 




385/311 


386/311 


326/311 


455 / 372 


374 / 330 


349/330 


410/347 


385/347 


479/347 


CXB MIN+NXB 3% RED 


fcr(keV) 


7.16 


6.61 


7.07 


7.09 


6.62 


5.95 


6.17 


2.53 


2.25 


z 


0.29 


0.33 


0.26 


0.19 


0.22 


0.27 


0.2 (FIX) 


0.2 (FIX) 


0.2 (FIX) 


norm 


650.2 


216.4 


69.4 


47.0 


28.6 


17.5 


5.8 


3.5 


1.4 


0^4-lOkeV 


94.93 


30.67 


10.65 


6.66 


3.91 


2.29 


0.78 


0.24 


0.09 


A-Vd.o.f 


387/311 


395/311 


323/311 


454/372 


364/330 


349/330 


402/347 


410/347 


450/347 


CONTAMI 10% ADD 


/tr(keV) 


6.62 


6.27 


6.75 


6.76 


6.13 


5.46 


5.00 


1.81 


1.31 


z 


0.29 


0.32 


0.25 


0.18 


0.22 


0.27 


0.2 (FIX) 


0.2 (nx) 


0.2 (FIX) 


norm 


650.2 


214.6 


66.2 


46.3 


27.0 


15.5 


5.2 


3.0 


1.0 


■^OJ-lOkeV 


91.87 


29.60 


10.10 


6.50 


3.61 


2.01 


0.62 


0.14 


0.03 


;t'Vd.o.f 


413/311 


365/311 


322/311 


455 / 372 


362/330 


347 / 330 


402 / 347 


402 / 347 


450 / 347 


CONTAMI 10% RED 


yfeT(keV) 


7.75 


6.98 


7.17 


7.58 


6.94 


6.13 


5.41 


1.93 


1.72 


Z 


0.30 


0.33 


0.26 


0.20 


0.23 


0.29 


0.2 (FIX) 


0.2 (FIX) 


0.2 (FIX) 


norm 


606.0 


201.8 


64.9 


44.3 


25.4 


14.7 


5.1 


2.9 


1.0 


'^0„4-10keV 


91.04 


29.46 


10.09 


6.49 


3.60 


2.02 


0.63 


0.14 


0.04 


;tVd.o.f 


536/311 


521/311 


326/311 


460/372 


392/330 


347 / 330 


406/347 


388/347 


450/347 


Stray light 


fcr(keV) 


7.06 


7.13 


7.26 


7.38 


6.99 


5.83 


4.09 


0.87 


0.87 


z 


0.29 


0.35 


0.28 


0.19 


0.23 


0.28 


0.2 (FIX) 


0.2 (FIX) 


0.2 (FIX) 


norm 


631.2 


210.9 


66.8 


46.8 


27.5 


16.7 


4.5 


2.3 


0.9 


04-lOkeV 


91.13 


30.75 


10.27 


6.58 


3.77 


2.08 


0.44 


0.03 


0.01 


//d.0.f 


383/311 


526/311 


325/311 


403 / 330 


367 / 330 


344 / 330 


406/347 


399/347 


494/347 



*:NormaUzation of the apec component scaled with a factor SOURCE-RATIO-REG/n,, from Table 2, 

Norm= SOURCE-RATIO-REG J n_.njjdV/(4w(l + z^)D\) x 10"^° cm"' arcmin"^, where Da is the angular diameter distance to the source. 

||: Photon flux in units of 10"* photons cm"^ s"' arcmin"^. Energy band is 0.4 - 10.0 keV. 

Surface brighmess of the apec component scaled with a factor SOURCE-RATIO-REG£2c from Table 2. 



ROSAT PSPC result as the input X-ray image (Knopp et al. 
1996). Table 4 shows the resultant photon fractions for each 
annular region. We can see that most (> 70%) of the observed 
counts come from the "on-source" region and 10-20% come 
from the adjacent region which is located on the cluster center 
side. Contribution from the bright central region (< 12.6') is 
less than 5% for the annular regions outside of 16.8'. 

4. Results 

4.1. Spectral fitting 

For spectral fitting, we estimate the efl'ective area for each 
annular using xissimarfgen (Ishisaki et al. 2007) assuming ayS- 
model (y8=0.54, rc=2.'97:Knopp et al. 1996) surface brightness 
profiles for the ICM of A3667. We carried out spectral fitting 
for the individual annular regions separately. 

We modeled the spectrum in each annulus as the sum of the 
ICM and the sky background components. Intensities, photon 
indices and temperatures of the CXB and the Galactic emission 
were all fixed at the values from table 3. We employed a single 



temperature thermal model {wabs x apec) for the ICM emis- 
sion of A3667. The estimated NXB spectrum was subtracted 
from each data set before the spectral fit. The interstellar ab- 
sorption was kept fixed using the 21cm measurement of the 
hydrogen column, A^h = 4.7 X lO^'^cm"'^ (Dickey & Lockman 
1990). The solar abundance in our analysis was defined by 
Anders & Grevesse (1989). We used XSPEC verl2.4.0 for the 
spectral fit. 

In the central regions within 25.'2, the free parameters were 
temperature kT, normalization norm and metal abundance Z of 
the ICM component. In the outer regions (> 25 .'2), we fixed 
the ICM metal abundance to 0.2, which is the typical value 
observed in cluster outskirts (Fujita et al. 2008). We used the 
energy ranges 0.35-8 keV for BI and 0.5-10 keV for Fl, respec- 
tively. The relative normalization between the two sensors was 
a free parameter in this fit to compensate for cross-cahbration 
errors. The spectra and the best-fit models for aU the annular 
regions are shown in Fig 3. The parameters and the resultant 

values are hsted in Table 5. 

We also examined the effect of systematic errors on our spec- 
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Fig. 4. (a) Suzaku BI and FI spectram of the 29.'4 - 33.'6 region fitted with the cie model in SPEX. To make it easy to see the "1 keV bump", we show 
the 0.35-2.0 keV spectrum. The lower panels show the residuals of a cie,2cie,pdem model fit to spectrum, (b) Emission Measure distribution. Each color 
shows the results of a different DEM model. All models show a twin EM peak around 0.9 keV and 4 eV, and the EM weighted temperature is 2. 1 keV, 
agrees with 1 kT model of the Suzaku (Mgenta: 1.92 keV, 29.'4 - 33.'6) and XMM (Gray: 1.9 keV for region 10, Finoguenov et al. 2010) results. 



Table 6. Fit results for the 29.'4 - 33.'6 region with 5 different thermal models. 



Method temperature (keV) norm* /d.o.f 

Ikf 1.92 ±0.38 1.23±0.11 388/347 

2kT 0.89+0.09,4.42+2.7 0.69+0.27, 0.57+0.33 360/345 

pdem - - 368/334 

c/em*(poly:n=8) - - 96/89 

dem*i8T) - - 114/81 



*: The dem can be used only as an independent model, 
therefore we subtracted all background components. 



tral parameters. We considered the systematic error for the 
NXB intensity to be +3% (Tawa et al. 2008), the fluctuation 
of the CXB is shown as CXBMAX and CXBMIN, and the 
systematic error in the OBF contamination is considered to be 
+ 10% CONTAMl+10% and CONTAMl-10% (Koyama et al. 
2007) in Table 2, respectively. 

To examine the influence of stray light on our spectral re- 
sults, we produced mock stray light spectra using xissim and 
fitted the annular spectra after this stray light spectrum was 
subtracted. To simulate the stray light spectrum by xissim, 
the spectral shape and the surface brightness distribution of the 
ICM are the input parameters. We assumed the uniform spec- 
trum that was observed in the central region (0'-4.'2: kT ~1.Q 
keV) and the y6-model surface brightness distribution described 
above. The resultant stray light corrected best fit values are 
shown in in table 5 in "stray light" row. We confirmed that 
the parameters after correction for the stray light effect stayed 
within the statistical and systematic errors of our spectral fit. 

4.2. Spectral features near the radio relic 

In the edge region of the radio relic (29.'4 - 33.'4), the sin- 
gle temperature model for the ICM does not give a good fit 
in the low energy range. As shown in Fig. 3(h), there is a 
residual feature around 0.9 keV which was already reported 
from Suzaku and XMM-Newton observations (Nakazawa et al. 
2009; Finoguenov et al. 2010). In this section, we look into the 
possibility of a multi-temperature plasma. 

To test the multi-temperature structure, we first included an 



additional thermal component {2kT model) to the model. The 
metal abundances were fixed at 0.2 times solar for both ther- 
mal models. We then obtained an almost acceptable fit with 

Id.o.f. - 360/345. The resultant high and low temperatures 
are /feThigh = 4.42 + 2.7 keV and kT]^^ =0.89 + 0.09 ke V, respec- 
tively. The single temperature fit gave kT - 1 .92 ke V, which 
was in the middle of these two temperatures. 

Additionally, we adopted differential emission measure 
(DEM) models using the SPEX spectral fitting pack- 
age (Kaastra et al. 1996). Several DEM models can account 
for complex emission measure (EM) distributions, and here we 
tried 3 models. The first DEM model we employed is a poly- 
nomial differential emission measure pDEM model. The sec- 
ond model is a polynomial method as dem (poly), which has 
been applied to several sources (Lemen et al. 1989; Stern et al. 
1995). In this method, the logarithmic emission measure is de- 
scribed as the sum of n Chebyshev polynomials as a function 
of logarithmic temperature. The last one is a Multi-temperature 
method as dem (multi T) , which is a sum of gaussian temper- 
ature distributions ; 

£M,o„i = y -^exp[-(log(r) - log(TDf/2cr^] (1) 
i y2na-i 

where EM, is the total emission measure of the component /, 
and Ti is the central temperature. 

Fig 4(a) shows the result of spectral fits with the pDEM 
model. The lower panel shows the residuals by fitting with 
\kT,2kT, and pDEM models. In the 2kT and pDEM cases. 
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Fig. 5. Radial profiles of (a) ICM temperature, (b) Surface brightness, (c)Abundance, (d)Deprojected electron density. The Suzaku 
best fit values are shown by black diamonds. In temperature profiles, gray diamonds show the results of XMM-Newton (Briel et al. 
2004; Finoguenov et al. 2010). Black dashed vertical lines show the approximate radial boundaries of the northwest radio relic. The 
Green and red dash lines show typical systematic changes of the best fit values due to changes of OBF contaminants and the NXB 
level. In temperature profiles, the crosses show an average profile given by Pratt et al. (2007) for A3667. In the surface brightness and 
deprojected electron density profiles, the black solid curve shows the /J-model obtained by ROSAT PSPC (Knopp et al. 1996). 



the residuals around 0.9 keV are less significant than for the 
IkT case. Fig 4(b) shows the resultant EM distribution based 
on the DEM analysis. In this analysis, we fixed the metal abun- 
dance to 0.2 solar (with the solar abundance table: Anders & 
Grevesse 1989) and the energy range from 0.5 to 5.0 keV. The 
result shown in Table 6. The resultant emission measure distri- 
bution indicates 2 peaks around 0.9 and 4 keV, which is very 
similar to the result of a 2kT fit. 

The emission measure weighted temperature becomes 2.1 
keV, which is consistent with the IkT results obtained with 
Suzaku and XMM-Newton (Finoguenov et al. 2010). In 
Finoguenov et al. (2010), most of the outer regions (No. 1, 
14 in Table 1) showed kT ~ I keV consistent with the results 
of the present analysis. Nakazawa et al. (2009) also reported 
a 0.9 keV component at the radio relic region using Suzaku 
data. Those results suggest either that (1) we observed the fore- 
ground Galactic component or that (2) the 29.'4 - 33.'6 region 
is a mixture of low and high temperature ICM. 



As discussted by Kaastra et al. (2004), a multi-phase gas 
makes an average temperature spectrum which is difficult to 
be distinguished by current X-ray satellites. To solve this type 
of problem by line diagnosis, a future X-ray satellite with a 
large grasp (FOVx Effective area) and good energy resolution 
is highly desirable. 

4.3. ICM properteis 

Recent studies from Suzaku have revealed the temperature 
structure in the outskirts of relaxed clusters, and suggested that 
ion-electron temperature equilibrium was not reached there 
(George et al. 2008; Reiprich et al. 2009; Hoshino et al. 2010). 
In the merging cluster Abell 2142, the temperature profile sug- 
gests some excess compared with relaxed clusters (Akamatsu 
et al. 2011). We have relatively little information about 
the ICM properties in the outer regions of merging clusters. 
However, our results on A3667 give us new information. Here, 
we use the ICM parameters based on our specti^al fits shown in 



8 



H. Akamatsu et al. 



[Vol. 



tables. 

Figure 5(a) shows the radial profile of the ICM temperature 
in A3667. The temperature within 15' (960 kpc) from the clus- 
ter center shows a fairly constant value of ~ 7 keV, and it gradu- 
ally decreases toward the outer region down to ~ 5 keV around 
the radio rehc (~ 21'). As mentioned in Section 1, Finoguenov 
et al. (2010) reported a significant drop of temperature and sur- 
face brightness in the radio relic region. Our temperature pro- 
file also shows a remarkable drop by a factor of 2-3 in this 
region even considering the systematic errors. 

We compare our temperature profile with the previous 
XMM-Newton result for A3667 (Briel et al. 2004; Finoguenov 
et al. 2010) and an average profile given by Pratt et al. (2007). 
The average profile was made from scaled temperatures for a 
sample of 15 clusters, and it is described by the following for- 
mula 

T/rx = 1.19-0.74/?//?200, (2) 

where Tx is the average temperature and 7?200 is the virial ra- 
dius. We used Tx = 7.0 keV and 7?200 = 34.' 1 for A3667. The 
temperature profile of A3667 shows good agreement with the 
previous XMM result for the radius range within the radio relic, 
as shown in Fig. 5(a). The observed temperatures of A3667 are 
somewhat higher than the average levels between 10' and 26'. 
Around the radio rehc region, the average profile predicts 4.0 
keV but the observed temperature is 5.3 keV, which is in excess 
by about 30%. 

The surface brightness profile shows a good agreement with 
the j8-model (fi = 0.54) up to the radio relic region. As shown 
in Fig. 5(b), the surface brightness also shows a drop by an 
order of magnitude across the relic. Figure 5(c) shows the 
metal abundance profile of A3667, which shows a drop from 
0.3 to 0.2 solar. Such a drop is seen in other merging or non- 
merging clusters based on previous studies by XMM-Newton 
(Matsushita 201 1; Lovisari et al. 2009). The abundance profile 
shows a significant excess around 0.1r2oo possibly due to the 
effect of the cold front. 

We calculated the deprojected electron density profile, based 
on the surface brightness results. The resultant profile is shown 
in Fig. 5(d) along with the previous XMM results (Briel et al. 
2004; Finoguenov et al. 2010), and the /?-model profile from 
ROSATPSPC (J3^QM,rc^2'.91). Our electron density profile 
and the XMM-Newton results are both in good agreement with 
the )8-model profile up to the radius of the radio rehc. 

Although the surface brightness of A3667 shows no signifi- 
cant excess over the /? model, the temperature profile is flatter 
and higher than the average profile. A similar feature was al- 
ready reported for another merging cluster, A2142 (Akamatsu 
et al. 2011). In A2142, the temperature profile along the 
merger axis shows shghtly faster dechne compared with the av- 
erage temperature curve for other relaxed systems studied with 
Suzaku. In the A3667 case, it is opposite because the tempera- 
ture shows a slower decline than in other relaxed clusters. The 
ROSAT and XMM images show an elongated shape along the 
merger axis direction. The temperature excess and elongated 
morphology suggest that the ICM is heated by the merger event 
and the time elapsed since then is not enough for full relax- 
ation. However, the ICM properties are only measured in the 
direction of the merger axis with this accuracy. Comparison 



with the ICM properties in the perpendicular direction would 
be important. 

5. Discussion 

Suzaku performed 3 pointing observations of Abell 3667 
along its merger axis. The temperature, surface brightness, 
abundance, and deprojected electron density profiles were ob- 
tained up to 0.9 times the virial radius (r2oo = 2.26 Mpc). The 
temperature and electron density show a significant "drop" 
around the radio rehc region, whose feature suggest a shock 
front. We evaluate the cluster properties (shock wave and en- 
tropy) and discuss their implications below. 

5.7. Possibility of a shock at the radio relic region 

Recent X-ray studies with Chandra showed clear evidence 
of a shock in IE 0657-56 and A520 (Clowe et al. 2006; 
Markevitch et al. 2005). However, there is still little in- 
formation about the nature of shocks in clusters. Recently, 
Finoguenov et al. (2010) reported a sharp edge in the surface 
brightness at the outer edge of the NW radio relic in A3667. 
They interpreted the feature in terms of a shock, and derived 
the Mach number M ~ 2.44 ± 0.77 and ~ 1.68 ± 0.16, based on 
the jumps of temperature and electron density in the radio relic 
region, respectively. 

The present Suzaku data show steep jumps in temperature, 
surface brightness, and electron number density at the region 
of the radio rehc. To verify the presence and properties of a 
shock, we estimate the ICM pressure. Fig 6(b) shows the re- 
sultant radial profile of the gas pressure in A3667. Black and 
gray curves show the results from Suzaku and XMM-Newton, 
respectively. Since the temperature and the surface brightness 
both show a steep decline, with a somewhat less significant 
jump in the surface brightness, we see a clear pressure jump of 
a factor of > 3. The spatial coincidence in the regions of the 
pressure drop and the radio rehc strongly suggests that there is 
a shock front in this region. 

We estimate the Mach number based on the Suzaku data in 
the same way as Finoguenov et al. (2010). The Mach number 
can be obtained by applying the Rankine-Hugoniot jump con- 
dition, assuming the ratio of specific heats as 7 = 5/3, as — = 

Ti , T2 and Pi ,P2 are the pre-shock and post-shock temperatures 
and pressure. C = n2/«i gives the shock compression. Using 
the observed jumps in temperature, electron density and pres- 
sure, we obtain M = 2.48 ± 0.46, 1.52± 0.10 and 1.87 ± 0.27, 
respectively. Table 7 compares our Mach numbers with the 
XMM-Newton values by Finoguenov et al. (2010). Suzaku re- 
sults agree well with the previous XMM results of At ~ 2.4 
with slightly smaller errors. 

Merger shocks in clusters were reported in the Bullet cluster 
(M ~ 3.0) and in A520 (M ~ 2.1). Although fliese Mach num- 
bers and the one in A3667 are broadly similar, their locations 
are very diflferent as pointed out by Finoguenov et al. (2010). 
In A3667, the shock is located very far from the cluster center 
at ~ 30' corresponding to 2 Mpc, in contrast to the shocks in the 
Bullet cluster and A520. Using the gas density in the outer re- 
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Fig. 6. Radial profiles of pressure, entropy, entropy ratio to theoretical value, and equilibration time after a shock passage, (a) (b) Black and gray 
diamonds show results of Suzaku and XMM-Newton, respectively. Dashed vertical lines show approximate radial boundaries of the northwest radio 
reUc. In the entropy profile (b), the black solid line shows K oc r'-^ given by Voit et al. 2003. (c) Ratio of the measured entropy with reflrect to the 
theoretical value (Kocr^-^). The horizontal scale is normalized by the virial radius. Black diamonds show Suzaku results for A3667, gray diamonds show 
A2142 results (Akamatsu et al. 2011), and gray dashed diamonds are A1413 values (Hoshino et al. 2010), respectively, (d) Ion-electron equiUbration 
time scale. Horizontal scale is normaUzed by the radius of the radio relic. The soUd curve shows the time after a shock heating assuming a constant shock 
speed V = 1400 km s"^ Black diamonds show the electron-ion equiUbration time calculated by equation (3). Gray soUd and dashed diamonds show the 
ionization time for n^t = 3 x lO'^andl x 10" s cm"^, respectively. 



Table 7. Mach number estimates for the A3667 radio reUc shock 



Mach number* 


Finoguenov et al. 2010 


This work 


M (by rie) 


1.68±0.16 


1.52 ±0.10 


M (by kT) 


2.43±0.77 


2.48 ± 0.46 


M (by Pressure) 


2.09±0.47 


1.87±0.27 



* We use the observed value as post shock region:25.'2 ■ 
and pre shock region: 29 .'4 - 33 .'6. 



■ 29.'4 



gion («e ~ 10 cm"^), the pre-shock sound speed is ~ Vss = 700 
km s"' . Using the shock compression C and this pre-shock ve- 
locity, we estimate a shock speed v shock = C • = 1360 ± 120 
km s~\ which is consistent with the previous XMM-Newton 



result {Vshock = 1210 ± 220 km s"^). This shock speed is much 
lower than in other systems (Bullet: 4500 km s"', A520: 2300 
km s ') and close to the speed of the cold front in the cluster 
central region Vcf = 1430 ± 290 km s"i (VikhUnin et al. 2001). 
Recent hydrodynamic simulations of merging clusters predict 
that cluster merger events generate cold fronts in the central re- 
gions and arc-Uke shock waves toward the outskirts (Mathis et 
al. 2005; Akahori & Yoshikawa 2010). These predictions are 
in good agreement with our observational results. 

Another possibility for the origin of the shock front is an 
infall/accretion flow from large-scale structures. Miniati et al. 
(2001) predict strong accretion shocks with Mach numbers up 
to 1000. Such shock waves caused by structure formation are 
also predicted in the Coma cluster (Ensshn et al. 1998). The 
Coma cluster shows a giant radio halo and a relic located far 
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from the cluster center, and it exhibits a temperature jump 
across the radio edge (Brown & Rudnick 2011). Although 
many simulations predict accretion shocks, the observational 
evidence is still very Umited and further sensitive X-ray obser- 
vations will be very important. 

Based on these observational and theoretical arguments, the 
shock at the northwest radio rehc of A3 667 seems to be more 
in favor of being caused by a merger event. 

5.2. Entropy profile 

The entropy of the ICM is used as an indicator of the en- 
ergy acquired by the gas. The entropy is generated during 
the hierarchical assembly process. Numerical simulations in- 
dicate that self-similar growth of clusters commonly shows en- 
tropy profiles approximated by r'-' up to r2oo (Voit et al. 2003). 
Since entropy preserves a record of both the accretion history 
and the influence of non-gravitational processes on the prop- 
erties of its ICM, entropy is a good indicator of the cluster 
growth. Recent XMM (Pratt et al. 2010) and Chandra re- 
sults on the entropy profile showed the entropy slope within 
rsoo, which is approximately 0.5 r2oo, to be in agreement 
with the predictions from numerical simulations. Recently, 
Suzaku has extended the entropy measurements close to ^200 
for several clusters, and showed a flattening or even a de- 
crease at r > 0.5 r2oo (Hoshino et al. 2010; Kawaharada et al. 
2010; Akamatsu et al. 2011; Simionescu et al. 2011). 

Comparing with relaxed clusters, we looked into the en- 
tropy profile of A3667. The entropy of the ICM is given by 
K = kTrie^^^. The entropy profile is shown, along with the 
XMM-Newton result (gray), in Fig. 6. The overplotted black 
line shows K oc r' '. The observed entropy slope is consistent 
with the theoretical value between 0' - 10'. The slope becomes 
flatter at r > 20' and suddenly drops around the radio relic re- 
gion. This drop suggests the existence of non-thermal heating 
around the relic. 

To compare the entropy profiles with the simulated slope of 
1.1 and with other clusters, we calculated the ratio of the sim- 
ulated and observed profiles. The resultant profile is shown 
in Fig. 6(c) together with the ones for an other merger cluster 
(A2142:gray) and a relaxed cluster (A1413:dashed). There is 
a clear deviation from the simulation in the range r > 0.4r2oo- 
Such a feature has been also reported in other relaxed clusters 
such as A1795, PKS0745, A1689 (George et al. 2008; Bautz 
et al. 2009; Kawaharada et al. 2010). The deviation from the 
simulated entropy curve in the cluster outskirts is apparently 
common to many clusters including merging and non-merging 
systems. Some physical conditions to explain the origin of the 
deviation were proposed (Simionescu et al. 2011; Hoshino et 
al. 2010). However, the interpretation of the observed features 
is not settled yet. 

5.3. The possibility of non-equilibrium ionization. 

Theoretical studies about a non-equilibrium state and an 
electron-ion two-temperature structure of the ICM in merging 
galaxy clusters have been carried out (Takizawa 2005; Rudd & 
Nagai 2009; Akahori & Yoshikawa 2010). Recent Suzaku re- 
sults suggest non-equilibrium features characterized by a dif- 
ferent electron and ion temperature in the cluster outskirts 
(Wong & Sarazin 2009; Hoshino et al. 2010). In merger clus- 



ters, shock heating and compression raise the X-ray luminosity 
in the outer regions, giving us a good opportunity to look into 
the non-equihbrium features and gas heating processes. 

We examine the possibility of a deviation between ion and 
electron temperatures. In A3667, the existence of a shock is al- 
ready suggested in the radio relic area (Finoguenov et al. 2010). 
Assuming that the shock propagates from the cluster center to 
the outer region, the region just behind the shock may have a 
higher ion temperature than the electron one, because it had not 
enough time to equilibrate. 

According to Finoguenov et al. 2010, we calculate electron- 
ion equilibration time and ionization timescales. The electron- 
ion equilibration time is estimated from the electron density 
and temperature of ICM. The ionization timescales can be esti- 
mated from relative line intensities. Because the Unes are very 
weak in the cluster outer region, we adopt the timescale for 
equilibrium ionization as n^f = 3 x 10'^ - 1 x 10^^ s cm"^ given 
by Fujita et al. 2008 indicating that the outer region attains ion- 
ization equilibrium. 

In fig 6(d), solid diamonds show the electron-ion equilibra- 
tion times of A3667 based on the observed temperature and 
density profiles. The gray and dashed diamonds in fig 6(d) 
show the resultant values of ionization timescales for riet = 
3 X 10'^ s cm^ and 1 x lO'"' s cm-', respectively. The horizontal 
axis of fig 6(d) is normalized by the radio relic radius. The sohd 
curve shows the expected time after the shock heating assum- 
ing a constant shock speed Vsi,„ck - 1400 km s"' . The relaxation 
and ionization timescales, considering electron-ion Coulomb 
colUsions, are both longer than the expected elapsed time after 
the shock heating around r > 0.7rshock- This suggests that the 
region near the radio relic would not have reached ion-electron 
equihbrium. This can be confirmed with future spectroscopic 
observations which wiU be able to determine ion temperature 
through line width measurements. 

6. Summary 

We observed Abell 3667 with Suzaku XIS and detected the 
ICM emission near the virial radius r2oo (2.3 Mpc ~ 34'). We 
confirm the jumps in temperature and surface brightness across 
the northwest radio relic region, which is located far from the 
cluster center (2 Mpc). Using the Rankine-Hugoniot jump 
condition and observed ICM pressure value, we evaluated the 
Mach number as At = 1.87 ± 0.27, consistent with the previ- 
ous XMM-Newton result (Finoguenov et al. 2010). The main 
results on A3667 are as foUows; 

• The ICM temperature gradually decreases toward the 
outer region from about 7 keV to 5 keV and then shows a 
jump to 2 keV at the radio relic. The temperature profile 
inside the relic shows an excess over the mean values for 
other clusters observed with XMM-Newton. 

• There is a distinct structure of the emission measure dis- 
tribution around the radio relic, characterized by two 
peaks around 0.9 and 4.0 keV. 

• The electron density profile shows a good agreement 
with the j6-model (/? - 0.54) and also shows a significant 
jump across the radio relic. 

• The significant jump in temperature and density profile 
across the radio rehc support the presence of a shock. 
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The estimated Mach number of the shock is A1 = 1 .87 ± 
0.27 based on the pressure jump. 

• The entropy profile within about O.Srioo follows \ pre- 
dicted by the accretion shock heating model. The profile 
shows a significant jump around the radio reUc, suggest- 
ing the region is not thermally relaxed. 

• Based on the relaxation and ionization time needed after 
a shock heating, we discuss the possibility of T, being 
higher than Te around the radio relic region. Because of 
the difl'erent parameter dependence of electron-ion equi- 
Ubration time and non-equiUbrium ionization, there is 
the possibiUty that not non-equilibrium ionization state 
but Ti and Tg are difl'erent. 

These results show that A3667 is a dramatic merging clus- 
ter, with the outskirts still in the process of reaching thermal 
equilibrium and the electrons in the northwest radio reUc be- 
ing strongly heated by their merger shock. Because of the ex- 
istence of a large cold front and presence of the radio relic, 
A3667 is a promising target for future X-ray studies such as 
with ASTRO-H. 
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ported by a Grant-in-Aid for Japan Society for the Promotion 
of Science (JSPS) Fellows (22-1582) and the JAXA ASTRO- 
H science team. H. A. deeply acknowledges hospitaUty in the 
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Appendix 1. Solar Wind Charge eXchange 

To check the contamination efl'ects of Solar Wind Charge 
eXchange(SWCX), we check the solar-wind proton flux during 
our Suzaku Observation Fig. 7. These plots are created by mul- 
tiplying the proton speed and proton density of public WIND- 
SWE data (http://web.mit.edu/space/www/wind.html). In the 
observation periods, proton flux has many flares. Because the 
ICM emission is very strong in the center region of the clus- 
ter, we can ignore efifects of SWCX. In the cluster outskirts, 
ICM emission is very weak, effects of SWCX is important. In 
this analysis, the light curve of the most outer region (OFFSET 
2) does not show time variation, therefore we assume efi'ect of 
SWCX was very smaU. 
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Fig. 7. (a) Proton flux (speed times density) measured by tlie WIND-SWE wiien SuzaJcu observed A3667. (b)-(d) Light curve of tlie Suzaku XISl 0.5 - 
2.0 keV band during each observations (OFFSET2, OFFSET 1 and CENTER). 



